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Abstract 

We study the Standard Model short-distance prediction for the mass and lifetime differences 
between the two neutral D meson mass eigenstates. We find that, despite a s /47r suppression, 
next-to-leading order (NLO) short-distance QCD corrections exceed the corresponding leading 
order (LO) amplitudes. For the lifetime difference, this stems from the lifting of helicity suppression 
of a light-quark intermediate state. We find y-Q is given by 2/nlo to a reasonable approximation 
but xd is greatly affected by destructive interference between xlo and £nlo- The net effect is to 
render yj) ~ id ~ 6 • 10~ 7 . Our NLO short-distance results, still smaller than most long-distance 
estimates, depend on the same two nonperturbative matrix elements of four-quark operators as in 
leading order. 
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I. INTRODUCTION 



i 



As such, there is 



Experimental efforts to detect D° — D° mixing are longstanding and remain an active area 
to this day. The theory of D° — D° mixing is relevant both in lending phenomeno- 

logical guidance to ongoing experimental work and in better understanding the workings of 
the Standard Model and of various New Physics scenarios . In this paper, we present 
new results — the perturbative QCD NLO contributions in the framework of the l/m c 
expansion for AI^ and AM^. The complex of D-meson phenomena presents a nontrivial 
theoretical laboratory for studying applicability of heavy quark methods. One can argue 
that m c 3> Aq CD justifies the use of heavy quark methods. However, the scale \i ~ M D lies 
in the meson resonance region, so QCD dynamics is clearly present 
inherently a degree of interest in the numerical aspect of our findings. Our calculation also 
touches on matters of principle, such as the degree of m q /m c suppression in ATp and AMp 
at NLO order. 

We begin by reviewing the theoretical context of D° — D° mixing. The mixing arises 
from AC = 2 interactions that generate off-diagonal terms in the mass matrix for D° and 
D° mesons. The expansion of the off-diagonal terms in the neutral D mass matrix to second 
order in perturbation theory is 

M '2 T ) 12 -2M^ {Dln - \ D) + 2M~ D \ Md — E n + ie ' (1) 

where 7i^ c=2 is the effective AC = 2 hamiltonian and Ti,^ c=1 is 

^ C=1 = % E Kfa PiWQi + c 2 Gu)Q 2 ] ■ (2) 

V 2 q,q> 

In H^ c=1 , the flavor sum on q, q' extends over the d, s quarks, 1 the quantities Ci^/i) are 
Wilson coefficients evaluated at energy scale /i, and are the four-quark operators 

Qi = (Vi c j)v-A ^4) V _ A and ^2 = (qiCi) v _ A (ujq' ) v _ A ■ (3) 

The first term in Eq. (JIJ represents AC = 2 contributions that are local at scale [i ~ M^, so 
it contributes to the M 12 (but not to the r 12 ) part of the mixing matrix. For example, in the 
Standard Model this term is generated by the contribution of the b quark. It can also receive 



1 In this paper, we work with m u — nid = 0. 
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a potentially large enhancement from new physics. The second term in Eq. comes from a 
double insertion of AC = 1 operators in the SM lagrangian, and it contributes to both Myi 
and T 12 . It is dominated by SM contributions even in the presence of new physics. At scale 
~ Mo, the contributions are from the strange and down quarks and these have relatively 
large CKM factors. By contrast the AC = 2 term is expected to give a negligible contribution 
(e.g. in the SM there is the severe CKM suppression \V ub V* b \ 2 /\V us V* s \ 2 = O(10" 6 )). Thus, 
we omit it henceforth. 

The off-diagonal mass-matrix terms induce mass eigenstates and D$ which are su- 
perpositions of the flavor eigenstates D° and D°, 

\D L , s )=p\D°)±q\D ) , (4) 

where \p\ 2 + \q\ 2 = 1. In the Standard Model CP violation in D mixing is negligible, as is 
CP violation in D decays both in the Standard Model and in most scenarios of new physics. 
We therefore assume in the rest of this paper that CP is a good symmetry, and adopt the 
phase convention [t| 

CV\D°) = -\D°) , (5) 

Then we have p = q, and |-Dl,s) become the CP eigenstates \D±) with CV\D±) = ±\D±). 
We then define the mass and width differences 

AM D = M D+ -M D _ and Ar D = T D+ - T D _ . (6) 

It is, however, customary to work directly with the dimensionless quantities, 

AM D Ar D 
= -= — , yv = -r=-, (7) 

Id 

where Td is the average width of the two neutral D meson mass eigenstates. 

The discussion thus far covers relevant background material. We conclude this section 
by addressing three particularly important additional points: 

I. Our calculation adopts an operator product expansion (OPE) Q, In the limit 
m c ^> A, where A is some soft QCD scale, the momentum flowing through the light 
degrees of freedom in the intermediate state is large. As such, an OPE is implemented 
by expanding the second term in Eq. (JTJ in series of matrix elements of local operators. 
For example, one writes for APd, 

Ar D = -2r 12 = -^l^(W\iJd'xT{H^ c =\x)H^\Q)}\D Q ) , (8) 



and expands the time ordered product in Eq. (JHJ) in local operators of increasing 
dimension (higher dimension operators being suppressed by powers of A/m c ). 

2. We calculate Ar^ by making direct use of work available in the literature ^(|, but 
not heretofore applied to D° — D° mixing. We then compute the mass difference AMrj 
from an unsubtracted dispersion relation, 2 

AM D (m 2 ) = -^pf°d S ^4 > ( 9 ) 
2n Jso s-m 2 c 

which follows from the analyticity of AI^ in the complex s-plane with a unitarity 
branch cut along the Tie s axis 

3. We expand all our expressions for x\ rm D and y^ in powers of the ratio z = m 2 s /m 2 c . 



II. ANALYSIS 



In what follows we compute LO and NLO contributions to y and then x, 

Vd = yho + Z/nlo and x D = x LO + x N lo • (10) 

We depict in Fig. ^ how QCD affects the D°-to-D° mixing amplitude: (a) the limit of no 
QCD corrections, (b) the LO component in which QCD dresses the interaction vertices, and 
(c) an example of a NLO correction. 




FIG. 1: D — ► D°: (a) No-QCD, (b) QCD-corrected vertices, (c) an example of NLO correction. 



The leading contribution to Ar^ in the 1 / m c expansion to D° — D° mixing comes from 
the dimension-six |AC| = 2 four-quark operators, 



Q = Ual^LCa U^^PlCp , Q S = U a P L C a UpPiCp 

Q' = Uofl^PLCp U^PlCo, , 0' s = u a P L c p u p P L c a 



(11) 



The tiny 6-quark contribution, neglected here, would contribute to a subtraction constant. 



where Pl = (1 + 7s)/2. One can use Fierz identities and equations of motion to eliminate 
Q' and Q' s in favor of Q and Qs- The resulting expression for AI"d is then 



F(z) (D°\Q^)\D°) + F s (z) (D \Q S W)\DP) , (12) 



12ttM d 

Coefficients F(z) and Fs(z) are defined as 

qq' 



f??(z) = f^\z) + ^fU™'(z), cm 



and similarly for Fg(z). Here £ 9 = V^* is a CKM factor for the intermediate s, d quarks, 
the {F^ qq (z)} functions are given in the discussion to follow and the {F^ qq (z)} are con- 
sidered later in our NLO analysis. As usual, the D°-to-D° matrix elements of Q and Qs are 
parameterized in terms of B-factors, 

(^IQI^ ) = lf 2 D M 2 D B D and (D°\Q S \D ) = -p 2 D M 2 D B^ , (14) 

where Ujf = B^ ] M 2 D /m 2 c . There are limits on the precision of Bd and B^ because the 
calculable short distance component most likely gives a negligibly small contribution. The 
most recent result for the quenched lattice calculation of £>d is reported in Ref. jl^ . 

A. Leading Order (LO) Contributions 

At leading order in a s , one finds for the ss intermediate state contributions to F(z) and 

F$ )ss (z) = 3v^4i (1 - z) Ff£\z) = 3VmIz~ (1 + 2z) 

F$ ss (z) = 2v^4i (1 - z) F s % ss (z) = 2v^4i (1 + 2z) (15) 

Fg> s \z) = 1(1 - Azf' 2 F s % ss (z) = -VT^Tz (1 + 2z) , 
and for the ds and sd contributions, 

Fi? ds (z) = 3(1 - zf (l + ^ F<g?{z) = 3(1 - z) 2 (1 + 2z) 



F^ ds (z) = 2(1 - z) 2 (l + |J F s % ds (z) = 2(1 - z) 2 (1 + 2z) (16) 
FT (z) = \{l - zf F s % d \z) = -(1 - zf (1 + 2z) . 



TABLE I: LO Values 





Without QCD 


With QCD (LO in z) 


With QCD (Exact) 


VLO 


-(2.9 -> 4.8) • 10~ 8 


-(5.6 -> 9.4) • 10~ 8 


-(5.7 -» 9.5) • 10~ 8 


XLO 


-(0.53 -» 1.05) • 10~ 6 


-(1.3 -> 2.3) • 10~ 6 


-(1.4 -» 2.4) • 10~ 6 



In addition, we have F { 



(0)dd, 



(0)ss 



(0). Insertion of Eqs. ffTT ^l . ([T5 |l . into Eq. © 



results in the following expression for the leading 0(z 3 ) contribution, 



Vho 



3vrr 



5, 



-5 



(17) 



where r^> ~ 1.6 • 10 12 GeV is the experimentally determined D° decay rate. The above 
expression for y^ Q agrees in the no-QCD limit of C\ — and C 2 — 1 with that found in the 
literature Since we expect 5B^ /2 > Bp, it follows that t/lo < 0. 

An expression for x^o is recovered by inserting AIYo into the dispersion relation of 
Eq. 0. One disperses in the variable m? c so that z = m? s /m 2 c — > m 2 s /s. The functions 
{F^\z)} of Eqs. (fI5 |l .(fIH |l are employed above the threshold for each intermediate state. 
Although the dispersion integral diverges separately for each of the ss, dd, ds, sd intermediate 
states, the flavor-summed expression for AM^ is rendered finite by GIM cancellations. All 
integrals are first evaluated analytically and the results are then expanded in powers of z. 
We find that the leading order in the ^-expansion for xlo occurs at 0(z 2 ), 

12 ,^2 f2 



X 



LO 



G 2 F mif 2 D M D 

Ss Z 



3vr 2 r 



D 



ri2n ^* ( r* 2 or* r< qri2\ n( s ) 



(18) 



As with yLO, we again regain the standard no-QCD result 0, |l3j]. Terms occurring at 



next-to-leading order in the z-expansion are straightforward to determine, and we find 



(z 4 ) 



(z 3 ) 
x LO 



G 2 F mif 2 D Mn, 2 4 



3vrr 



D 



B D (C 2 2 - AdC 2 - 6C 2 ) - ^-B^ (C 2 2 - 2C X C 2 - 3C 2 



GWJIM D _ 2 3 

Ss z 



3vr 2 r 



D 



B D (cf + 2C a C 2 + 3C 2 



h\z B 



'D 



25^(5) 



12 J 



5 



D 



c\ 



2C\C 2 — 3C 2 



Notice that at order £lo , there is now dependence also on hiz 
contributions are quite small relative to those of Eqs. (fTTj) . (fT%|) . 



(19) 
(20) 

-5. However, these 
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Numerical evaluations for y^o and xlo appear in Table HI The initial two columns display 
the leading z-dependences first with QCD turned off (cf Fig. E^a)) and then with QCD 
included via the Wilson coefficients of Eq. (J2J) (c/Fig. IHb)). The final column exhibits 
the exact LO results. The spread of values reflects uncertainties in input parameters (in 
particular, we have allowed for the range B^' / Bd = 0.8 — > 1.2. 

The collection of LO results in Table U gives rise to several interesting questions, but 
the most obvious one involves the tiny magnitudes. The main suppression arises from the 
presence of z 3 ~ 2 • 10~ 7 in y^o and z 2 ~ 4 • 10 -5 in xlo, even though the expansions 
for Fij(z) and Fsij(z) begin at 0(1). Such 0(1) contributions would be enormous, but 
they are in fact cancelled away as are O(z) terms. As a result, ?/lo an d x^o are rendered 
tiny. A numerical by-product of the dependence j/lo ~ 0(z 3 ) and x^o ~ 0(z 2 ) is that 
I^loI 3> |2/lo I - There is of course a corresponding physics explanation. In the diagrams of 
Fig. HI the 6-quark contribution is severly CKM suppressed, so only the light s, d quarks 
propagate on internal legs. Since the mixing amplitude will vanish in the = m s = limit, 
the breaking of chiral symmetry and of SU(3) flavor symmetry play crucial roles. Thus, a 
factor of m 2 comes from an 577(3) violating mass insertion on each internal quark line and 
another from an additional mass insertion on each line to compensate the chirality flip from 
the first insertion. This mechanism of chiral suppression accounts for the z 2 dependence of 
xlo m addition, y LO requires yet another factor of m 2 oc z to lift the helicity suppression 
for the decay of a scalar meson into a massless fermion pair. 



B. Next-to-Leading Order (NLO) Contributions 

Any way of reducing the chiral and helicity supression in x and y should lead to an 
enhancement. In principle, there are both nonperturbative and perturbative ways to achieve 
this. 

One might associate nonperturbative effects with the presence of quark condensates in the 
QCD vacuum 0,1^. These contributions (suppressed by powers of l/m c ) lead to chirality 
flip the same way mass insertions do, but have an intrinsic scale of A ~ 1 GeV ^> m s . In the 
realistic case of not-so-large m c , such power suppressions are not always sufficient to ensure 
the smallness of higher order contributions. Therefore, Eqs. (117)1 and ()18j) cannot contribute 
to leading order in the dual expansion in m s and l/m c if higher order terms in the l/m c 
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expansion contain lower powers of z than do xlo and y-^o- It has already been shown that 
this is indeed the case 0,0]. 

There are also perturbative QCD corrections to x and y, but these have heretofore not 
been given serious consideration due to the negligibly small LO values for the D° — D° 
mixing parameters. Even taking into account large scale dependence, the LO result gives a 
tiny contribution. This has stimulated a shift of attention towards the computation of the 



15] 



long-distance sector with varying degrees of model dependence jljj, |1J, 

But due to their milder dependence on m s , the higher order QCD corrections might be 
able to give relatively large contributions [ifj]. This occurs, for example, in the c — > 57 short 



distance amplitude, which receives a huge QCD correction |17| . In this paper, we consider 
a specific means by which the helicity supression in y can be lifted — a perturbative gluon 
correction (e.g. as in Fig. dc)). Having a perturbative gluon traversing the graph for the 
correlation function is the same as a well-known effect of lifting of helicity supression which 
follows from having three particles in the intermediate state instead of two 3 . The addition 
of the 'intermediate-state' gluon can lift one power of z, which characterizes the helicity 
suppression in y. Also, the relative lightness of m c implies that higher order perturbative 
QCD corrections are suppressed by a relatively large factor of a s (m c ) ~ 0.4. It is therefore 
expected that the NLO corrections to y should dominate the LO result. Moreover, the 
existence of a dispersion relation implies that x might well be enhanced at NLO. 

In order to systematically include the effects of intermediate-state gluons, a complete 
calculation of NLO corrections to D° — D° mixing is needed. The NLO corrections to lifetime 
difference y can be readily computed. All the relevant NLO contributions to F(z) and Fs(z) 
for two massive quarks and one massive, one massless quark can be found by adopting 
the formulas in Refs. [10] (which considered the case of B s — B s mixing) to computing 
F^p 9q of Eq. ()13j) . That calculation has been performed in the NDR-scheme (dimensional 



regularization with anti-commuting 75 and MS subtraction). We shall not present explicit 
formulas for the {F-p 9q (z)} and {F^} 99 (z)} functions as they are rather cumbersome. 
Scale dependent quantities used in our numerical work and evaluated at /1 = 1.3 GeV 



3 This mechanism leads to the prediction that the rate for the weak radiative decay B —> ^ev is much larger 
than the rate of weak leptonic decay B — > ev 
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TABLE II: NLO Values 



(z 2 ) 
^NLO 


(z s ) 


2/NLO 


^NLO 


(PENG) 

"nlo 


(2.2 -► 6.3) • 1(T 7 


(1.7 2.8) • 1(T 7 


(3.9 -► 9.1) • 1(T 7 


(1.7 3.0) • 1(T 6 


(0.6 -► 0.8) • 10~ 9 



were: 



m c = 1.3 GeV , B D = 0.82 , d = -0.411 , C 2 = 1.208 , a s 



0.406 



(21) 



The value for B D at scale fj, — 1.3 GeV is obtained by referring the lattice determination at 
/i = 2 GeV and employing the scale invariant quantity S D , 

a s (/i ) 



S D = 5 D (/zo)[a s (/i )r 6/25 



1 + 



47T 



with J 4 ~ 1.792. Also we allow for a range of the ratio B D /B D . 



Ja 



(22) 



Using the {F^ qq '(z)} and {Fgg 9 (z)} functions, we have calculated 2/nlo exactly and 
also have expressed it in terms of a power series in z. The leading term is 0(z 2 ), 



( 2 ) G 2 F mlfpM D 2 m s , 
Z/nlo — 



37rr D 
2 D 



47T 

'8tt 2 



25\ 



77 
6 



C 2 2 + 14 CiC 2 + 8C 
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C 2 + 20 dC 2 + 32 C 2 



(23) 



and the corresponding 0(^ 3 ) contribution is 
1G\m\p D M D - a s 



(3) 

Z/nlo — 



x LBr 



3vrr 



4tt 



'77 103 

l5 + 71nz) C 2 - ( — + hiz) CiC 2 - (18 + 58 In z) C\ 

9 3 



-5 



o 



28 
3 



^ + 61n^) C\ + ( ^ 



|_H? ln ,) ClC2 _(^ + 581nz) ^ 



(24) 



The numerical results, displayed in Table II, reveal that t/nlo is almost an order of magnitude 

(3) (2) 

larger than y^o and that the subleading term y^Lo ^ s smaller than y^LO but not at all 
negligible. 

The corresponding expression for xnlo has, as before, been obtained by means of a 
dispersion relation. We evaluated the dispersion integral numerically to obtain the value 
presented in Table II. As regards an analytical expression for xnlo, the intent was again 
to by first exactly perform the dispersion integrals and then expand each contribution in 
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a z power series. It turned out possible to do this for the dd, ds and sd intermediate 
states, but not for ss. It is, however, nonetheless useful to have an approximate analytic 
representation for Xnlo- By exploring a variety of approximation techniques, we found 
the expected 0(z 2 , z 2 \nz) leading behavior for £nlo but encountered scatter in the 0(z 2 ) 
coefficients, although less so for the 0(z 2 In z) coefficients. Upon accepting the latter and 
fitting the 0(z 2 ) coefficients to the numerical evaluations of individual dispersion integrals, 
we arrived at the 'effective' formula: 

G 2 F m 2 J 2 D M D 2 a s 2 
XNLO " 3^f^4^ 

x(b d [(11.3 -4.1 In z) Cl + (49.2 + 15.8 In z) C X C 2 + (37.9 + 10.7 In z) C 2 ] (25) 

--Bl [(37.9 + 2.2 In z) C 2 + (-33. + 81.8 In 2) C X C 2 + (32.0 + 125.3 In z) C 2 

8 l 

This relation, although approximate, is nonetheless useful in understanding the magnitude 
of the various contributions to £nlo 

Since the NLO results found for the box contributions are larger than their LO counter- 

(p) 

parts, we consider here for the sake of completeness the NLO penguin contribution ?/n LO to 
the width difference. We have 

i/ffio = - 4G \ m jg MD e £ a (bu + ae) + . . . . (26) 

The result shown in Table m clearly shows the penguin amplitude for t/^lo ^ s negligible 
compared to the box contribution. The mass splitting x^ Q is likewise 0(z 3 ) and hence 
negligible. 



III. CONCLUDING COMMENTS 



We have calculated LO and NLO contributions to the leading dimension-six component 
in the OPE for D° — D° mixing. Numerical results appear in Table IJfor LO and Table ITT1 for 
NLO. As a partial check of our analysis, we found our results (in cases of overlap) to agree 
with work carried out previously. Our formulae for x and y involve not simply expansions in 
l/m c , but rather combined expansions in m s (rrid is negligible), a s , and l/m c . As a technical 
aside, we performed the calculations at scale m c ~ 1.3 GeV. 

The two most noteworthy numerical features found for x and y are: 
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1. They are small at LO and even at NLO. This is because z = m 2 /m 2 is small and the 
leading dependence on z is found to be 

2/lo ~ z 3 x LO ~ z 2 y NLO ~ z 2 x NLO ~ z 2 . (27) 

Although contributions from individual intermediate states are not small, CKM factors 
cancel away the 0(1) and O(z) components. 

2. The NLO terms are larger than the LO terms. This requires somewhat more expla- 
nation, especially since NLO amplitudes contain the small perturbative QCD factor 
a s /47T. As regards the dimensionless width difference y&, the ratio of leading terms in 
the z expansion is 

>™g. = £i x - x * ~ 0.03 x 169 x (-0.73) ~ -4 . (28) 

JUS 47r Z wi } 

In the above W^ nlo ^ /W^ L °^ is the ratio of terms containing the Wilson coefficients in 
Eqs. (1171) , (|23j) and for definiteness we have considered the case = 0.8B D . Eq. 

(z 2 ) (z 3 )i 

shows that I^/nloI exceeds ly^o I because the extra factor of z overwhelms the a s /47r 
suppression. We have already discussed the physics of this - the helicity suppression 
mechanism which affects any LO qq intermediate state is removed via the presence of a 
virtual gluon in the NLO qqG intermediate state. Also, the difference in sign between 

and Vho arises from the factor W^ L0 yW^ L0 \ 
Since, to leading order in z, £nlo an d £lo both behave as z 2 , something else must 
account for the result |xnlo| > I^loI- From Eq. (fTHJ) and the approximate formula 
Eq. (|2*5j) . we have 

T^T - x - -0-03 >< <«-4) - -1-3 , (29) 

where W^ h °^ and are again the contributions from the Wilson constants and 

their coefficients. In this case, the suppression in a s /4n is overcome by the large size 
of Wj NL °) /W^ LO \ In particular, the largest contributor to W / "j NLO) is from the 
term in Eq. (|25|). roughly equally between log and non-log terms. 

3. We conclude that, citing just central values, the net effect of the short distance con- 
tributions is 

Vd = Vlo + Vnlo - 6 ■ 10~ 7 , x D = x LO + x NLO ^ 6 ■ 10~ 7 . (30) 
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In brief, yr> is given by t/nlo to a reasonable approximation but id is greatly affected 
by destructive interference between xlo and xnlo The net effect is to render t/d and 
xd of similar magnitudes, at least at this order of analysis. 

D° — D° mixing thus provides a concrete example of a well-defined observable for which 
NLO perturbative QCD corrections dominate the LO result. Will it follow that the NNLO 
contributions are larger still? Of course, one cannot know without doing the calculation. 
We feel, however, it may not necessarily be the case, at least as regards the width difference 
y. The three-particle qqG intermediate states were able to lift the helicity suppression 
experienced by qq intermediate states. In passing to the NNLO sector, there is no analogous 
suppression factor to be lifted. Of course, there is always the possibility that large numerical 
coefficients can overturn the a s /4n counting. 

The question remains - just how large is D° — D° mixing? Evidently, it is still not 
possible to provide a definitive theoretical answer and experiment will presumably decide 
the issue. On a relative basis, our 'short- distance' numerical results are smaller than most 
'long-distance' estimates (although the model dependence and uncertainty present in even 
modern and improved versions of the latter is less significant here). Experimentalists might 
find it useful to interpret our numerical NLO values as lower bounds to and xd- 

We conclude by considering our analysis in the context of operator product expansions. 
As we have seen above, the prediction of x and y is a result of expanding the correlation 
function Eq. (jHJ) in terms of three 'small' quantities, z, A/m c , and a s . Since the first quantity 
is significantly smaller than the other two, the structure of the series is rather different from 
other (usual) applications of the OPE, e.g. B° — £?° mixing or 6-hadron lifetimes [lsj ]. 

Working with this combined expansion, we computed the leading contribution originat- 
ing from matrix elements of dimension-six operators. These matrix elements are commonly 
parameterized in terms of the two nonperturbative parameters, Bp and Bp. The applica- 
tions of techniques of lattice and QCD sum rule evaluations of these operators can hopefully 
further improve the precision of our prediction. 

At higher orders in this expansion one would need to take into account 0(z 3 ^ 2 ) cor- 
rections (multiplied by about a dozen matrix elements of dimension-nine operators) and 
O(z) corrections (with more than twenty matrix elements of dimension-twelve operators). 
This would introduce a veritable multitude of unknown parameters whose matrix elements 
cannot be computed at this time. Simple dimensional analysis j^] suggests magnitudes 
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xd ~ Ud ~ 1CF 3 , but order-of-magnitude cancellations or enhancements are possible. How- 
ever, any effect of higher orders in l/m c or a s (m c ) which could render the result to be 
proportional to z n in the lowest possible power n = 1 [15] would presumably produce a 
dominant contribution to the prediction of x and y. 
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